Manipulating Higher Partial- Wave Atom- Atom Interaction by Strong 

Photoassociative Coupling 
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We show that it is possible to change not only s-wave but also higher partial wave atom-atom 
interactions in cold collision in the presence of relatively intense laser fields tuned near a photoas- 
sociative transition. 

PACS numbers: 34.80.Qb, 34.50.Cx, 32.80.Qk, 34.20.Cf 



Ability to control particle-particle interaction is im- 
portant for exploring quantum physics of many-particle 
systems in various interaction regimes. Ultracold atoms 
offer a unique opportunity for such explorations with un- 
precedented control over atom-atom interaction. There 
are two methods of manipulating interaction in cold 
atoms. The most popular one is magnetic field Fcshbach 
resonance (MFR) |l[ which has been extensively used to 
tune s-wave scattering length over a wide range. This has 
facilitated the recent demonstration of s-wave fermionic 
superfluidity in strongly interacting atomic gases Q. In 
fact, MFR has become an essential tool in experimen- 
tal investigations on the effects of large s-wave scattering 
length on the properties of atomic Fermi gases 0] and 
Bose-Einstein condensates (BEC) Q. The other method 
of modifying atomic interaction is optical Fcshbach reso- 
nance (OFR) proposed by Fedichev et al. [B| and imple- 
mented in recent experiments @, 0, H . While MFR relies 
on magnetic effects of Zeeman and hypcrfinc interactions, 
OFR uses off-resonant continuum-bound optical dipole 
transitions. In the case of resonance or near-resonance, 
OFR can lead to photoassociation (PA) Q of two atoms 
into an excited molecule. Recently, p-wave MFR [l(| 
in fermionic atoms has been observed. Enhanced scat- 
tering in higher partial waves by magnetic-field induced 



dissociation of Feshabch molecule has been shown 11]. 



There is a proposal 12j for generating anisotropic in- 
teraction by static electric field. Both the methods of 
magnetic and optical Feshbach resonances are so far pri- 
marily used to tune s-wave scattering length in ultra- 
cold atoms. To go beyond s-wave physics of cold atoms, 
it is now essential to devise methods of controlling p-, 
d- and other higher partial-wave interactions. This is 
particularly important for testing models of unconven- 
tional superconductivity or superfluidity in atomic Fermi 
gases. Superfluidity and superconductivity are related 
phenomena. Conventional low temperature supercon- 
ductivity can be explained by Bardeen-Cooper-Schrieffer 
theory which is based on s-wave Cooper-pairing. It is as- 
sumed that higher partial-wave interactions can lead to 
unconventional and high temperature superconductivity. 
Studies on Fermi superfluidity in cold atomic gases with 
controllable p- and d-wave interactions will help us to de- 



velop new insight about high temperature superconduc- 
tivity which requires a proper theoretical understanding. 

Here we show that it is possible to change not only 
£ = (s-wave) but also nonzero partial-wave scatter- 
ing amplitudes of two cold atoms by OFR with a rela- 
tively intense laser field. At low energy, the light-shift (or 
Stark-shift) due to laser-induced free-bound coupling can 
greatly exceed the spontaneous as well as stimulated line 
widths of excited molecular state. An intense PA laser 
can set in two photon processes in which one photon will 
cause PA transition from continuum to bound level and 
another photon will induce stimulated transition back to 
the continuum. If the light-shift largely exceeds the stim- 
ulated line width, then even when PA laser is tuned near 
the unperturbed (without Stark shift) excited molecular 
level, the formation of excited molecule becomes unlikely 
due to large light-shift. In such a situation, s-wave scat- 
tering wave function can be made to couple to p-wave 
or even d-wave scattering wave functions depending on 
the coupling of the molecular axis with the electronic or- 
bital and spin angular momentum. Furthermore, it is 
possible to enhance £ ^ partial wave scattering ampli- 
tudes with multiple strong laser fields causing continuum- 
bound PA coupling with appropriate rotational states of 
an excited vibrational level as illustrated in Fig.l. In 
this context, recent experimental results Q on the in- 
tense laser field PA of ytterbium may be of relevance. In 
molecular bound-bound spectroscopy, it is known that 
the rotational states of a molecule can be excited by in- 
tense laser fields [l]| , but exciting higher partial waves in 
continuum states by continuum-bound PA spectroscopy 
has not been considered so far. 

Let us consider that the scattering state of collision 
energy E — h 2 k 2 /(2/x) (where n is the reduced mass) of 
two colliding ground state atoms is coupled to an excited 
molecular state characterized by v vibrational and J rota- 
tional quantum numbers. The electronic orbital (L) and 
spin (S 1 ) angular momentum of the excited diatom are 
coupled to the diatomic axis according to either Hund's 
case (a) or (c). In Hund's case (a), A and £ which are 
the projections of L and S, respectively, on the internu- 
clear axis are two good quantum numbers and so is their 
sum = A + E. In Hund's case (c), the projection £1 of 
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FIG. 1: A schematic diagram of ground and excited poten- 
tials, scattering and bound states, rotational levels and rele- 
vant PA transitions for modifying £ > partial wave scatter- 
ing amplitudes. An intense laser (double arrow red line on 
the left) is tuned near J — 1 rotational state (of a particular 
vibrational level v) which can be accessed by PA transitions 
from s-wave and also at least the next nonzero (p-wave) scat- 
tering state. In the strong coupling dispersive regime with 
large light shift (see text), p-wave scattering amplitude will 
get modified due to its indirect coupling with s-wave scat- 
tering state. This modification can be probed by sending a 
weak probe laser (single arrow blue line) resonant with J = 3 
PA transition. The modification of p-wave can further be en- 
hanced by applying another intense laser (double arrow blue 
line on the right) tuned near J = 3 PA transition. 



the total electronic angular momentum J e = L + S is a 
good quantum number. The angular state of diatom can 

be written as | JflM) = i J \J^^'D^ I ) n (f) where M is 
the z-component of J in the space-fixed coordinate (lab- 
oratory) frame. T> (f) is the rotational matrix element 
with f representing the Euler angles for transformation 
from body-fixed to space-fixed frame. For Q = 0, T^Mrii^) 
reduces to spherical harmonic Yjm- The dressed state of 
a bound level (v, J) coupled to a continuum can be writ- 
ten as 

9e = ^ WO ! mM) ! e)eiec 



J dE'fa, 



£m e 0) | g)elec (1) 



where ^jmW is the radial part of the bound state, | 
a)eiec denotes internal electronic part of the excited (a = 
e) and ground (a = g) molecular states, ipE'im e (r) rep- 
resents energy-normalized partial wave scattering state 



with collision energy E' . Here |/?b'| 2 denotes density of 
states of the unperturbed continuum. The internal elec- 
tronic states | a) are the functions of electronic coordi- 
nates and have parametrical dependence on the internu- 
clear separation r. In electric dipole approximation, the 
interaction Hamiltonian is Hint — X)i=i 2 ^L^-di where 
di = —evi is the dipole moment of i-th atom whose va- 
lence electron's position is given by with respect to 
the center of mass of this atom. Here e represents an 
electron's charge, El is the laser field amplitude and tt 
is the polarization vector of the laser. In the absence 
of hyperfine interactions, the total Hamiltonian in the 
center-of-mass frame of the two atoms can be written 
as H = H elec (r 1 ,r 2 ;r A ,r B ) ~ f^V 2 - |^V| + H mt . 
where H e i ec includes terms which depend on electronic 
coordinates. Here ya and yb represent the position 
vectors of the nuclei of atoms A and B, respectively, 
V r and Vfl denote the Laplacian operators correspond- 
ing to the relative coordinate r = ya — ?b and the 
center-of-mass coordinate R = (r^ + yb)/2. From time- 
independent Schrodinger equation H^e = E^e, using 
Born-Oppenhcimer approximation we obtain the follow- 
ing coupled equations 



2fi dr 
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imp 



E 



/ J lUm e ;JM<PvJM 
M 



(3) 



where Bj(r) = ti 2 / (2(ir 2 )[J(J + 1) - D, 2 } is the ro- 
tational term of excited molecular bound state and 
Bi(r) — h 2 /(2/j,r 2 )£(£ + 1) is the centrifugal term in 
collision of two ground state atoms and ipE£m e (f) = 
J" , PE'tpE'ime (r)dE' . The free-bound coupling matrix 
element is A JM , e „ le = ( JMQ, | (0 e (ri,r 2 ;r) | H mt \ 
<f)g[ri, r 2 ; r)) | £m^S). The molecular electronic wave 
functions <^> a (ri,r 2 ; r) = (ri,r 2 ;r | a) can be constructed 
from the symmetrized (or antisymmetrized) product of 
atomic orbital of the two atoms using Movre-Pichlcr 
model [14] which also provides the long-range part of 
adiabatic potentials V a . We have here introduced a term 
hry/2 corresponding to the natural line width of the ex- 
cited molecular state in order to take into account the in- 
elastic process of natural decay of the bound state. V e (r) 
goes as — 1/r 3 for r — > 00 while V g (r) behaves as — 1/r 6 
in the asymptotic regime. The excited state potential 
V e supports several bound states. Here S = ll>l — ^>A is 
the frequency off-set between the laser frequency lol and 
atomic resonance frequency uja- The coupled Eqs. ([2|) 
and can be solved by the method of Green's function. 
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Let 4>vj be the bound state solution of the homogeneous 
part (A = 0) of © with rovibrational energy E v j. Note 
that we have here removed the subscript M in the la- 
belling of wavefunctions for simplicity. The correspond- 
ing Green's function can then be written as 



G v (r,r r ) = - 



ifvy/2 



AE, 



vj 



(4) 



where AE v j = HS+E—E v j. We can express the solution 
of equation © in the form 



4>vJM(r) 



dE'Y,Aj M -4m^E4lj{r) (5) 



where 



A 



-JM;imt — 



J dr'Aj M -£m e (r')(f>° j(r')ip Etme (r 1 ) 
AE vJ + ihj/2 



(6) 



The Green's function for the homogeneous part of © 
can be constructed from the scattering soultions. Let 
i^Ei eQ ( r ) an d ( r ) represent the regular and irregular 
scattering solutions in the absence of laser field. ^eY 9 ( r ) 
vanishes at r — while i^Ei[ r ( r ) is defined by boundary 
condition at r — * oo only. Asymptotically, they behave as 

69 ( r ) ~ 3l cos Vi - n i sin Vi an d ^e" ( r ) ~ n t cos % + 
ji smrji, where r\g is the background phase shift of i-th 
partial wave in the absence of light field and jt an d ne 
are spherical Bessel and Neumann functions. According 
to threshold laws, as k — ► we have r\i ~ k 2e+1 for £ < 
(n— 3) /2, otherwise rj£ ~ k n ~ e ; with n being the exponent 
of the inverse power-law potential at large separation. 
The Green's function [l5| for the scattering wave function 
can be written as 

Ke(ry) = -7r[^r i '(^r(r , ) + i^ refl (r)^(r')] 

r' >r (7) 

Wry) = -^[^rnr0^r(r) + z^; e9 (r)^f 9 (r')] 

r' <r (8) 

Substituting equation ([5]) into equation ([3]) and using 
/Q(r, r'), we have 

ipEe me (r) = exp(irji)ip% e + ^ A .JM-i> me ,{E) 

I'm.iM 



E 



shift 

Jl ~ 



dr'dr<l> v j{r')AjM;em e (r')Re[AC/(r', r)] 



x A emi; jM(r)4>°j(r). 



(11) 



where fjM;tm t 

Aj = E 



shift 



exp(%) / 4Pj(r')AjM;£m e ( r ')' t l , %e( r ')dr' , 
I <P vJ (r)AjM;tm e {r)i>%t{r)dr\ 

f.JM;tmf 

e,m ( ,M AEvj+ihy j2-Ef %st +ihT ., /2 ' 



and 
Here 



J - Z^.m^M^j/ and Tj = E£,m f ,M r^- 

Using the asymptotic boundary conditions of regular 

and irregular scattering wave functions, the scattering 

T-matrix in the presence of light field can now be written 

as 



T t = - [exp(2z % ) - 1] 
Li 



exp[2irn] 

~^2l'M ftme;JMfjM;£'m e , 



H5 + E- E vJ + ihj/2 - Ej 



shift 



ihTj/2 



(12) 



The partial wave S-matrix element can now be obtained 
from the relation Si = 1 + 2iTp. Only the second term 
on right hand side (RHS) of the above equation con- 
tains the effect of light field. This equation reveals that 
the T-matrix element for a partial wave £ ^ can be 
modified by its indirect coupling with £ = via the 
excited rotational state. The modification for £ = 1 
is mainly due to the two-photon transition amplitude 

tfi = ft'=o->J=ifj=i->t=i- Tnc shift °f Eq. dTTJ> in- 
volves the real part of the Green's function /Q(r, r') and 
the bound wave functions at two space points r and r' . 
Thus the shift depends on the radial correlation between 
continuum and bound states. In the limit k — > 0, the 
shift becomes independent of collision energy for all par- 
tial waves. For large k the shift will be vanishingly small. 

For numerical illustration, we consider a model system 
of two cold ground state Na atoms coupled to v = 48 vi- 
brational state of l g molecular potential by a laser field. 
Recently, several (up to J = 6) sharp rotational lines of 
this vibrational state have been observed in PA spectra 
with a strong laser field [161 ]. The outer turning point of 
the excited v = 48, J — 1 level lies inside the centrifugal 
barrier of I ^ of the ground continuum. Therefore, 
nonzero partial waves are not expected to contribute sig- 
nificantly to the PA transition amplitude at low temper- 
ature in the weak-coupling regime. This situation can be 
contrasted to the PA spectroscopy of higher rotational 
states where transitions occur outside the barrier region 
171 ] . The results of our numerical calculations as tabu- 



x f fCe(r,r')Ae me . t j M (r')<p v j{r')dr' (9) lated in Table-I show that the light shift Ef ift can ex- 
J ceed the stimulated line width Y i bv more than one order 



where ij) El = ^eY 9 • On substitution of equation §§§ into 
([6]) and after some algebra, we obtain 



.4 



J,M;l,m t 



fj,M;lm e 
1 



iKTj t /2)A 



AE vJ + ih-f/2' 



(10) 



ceed the stimulated line width Tj by more than one order 
of magnitude when the PA laser intensity is as high as 10 
kW/cm 2 . The natural line width 7 of the ro- vibrational 
states is of the order of 100 kHz [18[. However, the light 
shift remains much smaller than the rotational energy 
spacings Aj = E v j + \ — E v j. Since the background 
(without laser) phase-shift rje. — * as the collision en- 
ergy E — > 0, we can approximate tf^> ~ sf^j^ji 1 /2 as 
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TABLE I: Numerically calculated partial energy shirts 
and partial stimulated line width Y ji (in unit of MHz) for 
PA laser intensity I — 10 kW/cm 2 and collision energy E = 
50/iK. Also given are the rotational energy spacings Aj = 
E v j+i — E v j (in unit of GHz) for a few lowest J values. The 
total shift and stimulated line width for J — 1 are E^"'* = 
62.70 MHz, Ti = 2.87 MHz, respectively. 
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£^ l/ *(MHz) 
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J A j (GHz) 
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-14.22 
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1 1.56 
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-17.20 
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-15.28 


10" 4 


3 3.78 
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-16.00 


10 -8 


4 4.48 



real quantity unless the laser introduces a phase. It then 
follows that the elastic scattering will be predominant if 
the condition (AE vJ ~ Ef ift ) » h(-y + Tj) is fulfilled. 
As k — > 0, in the leading order in k the ratio of the laser- 
induced change in p-wave T-matrix element to that in 
s-wave one is given by \/T j = u=i/^ j=u=o — 0.28. We 
can write an energy-dependent €-wave scattering length 
as ae — — = a° + af, where a° is the background scat- 
tering length and af denotes the laser-modified part of 
ai. Note that the scattering length ai ( for I ^ 0) as 
defined here differs from the standard definition in scat- 
tering theory However, ag as defined here can be 
related to the standard £-w&ve scattering length by us- 
ing the behavior of ag in the limit k — > and thereby 
can be compared with the results of Ref. [2(j. When 
S vJ — hS + E — E vJ > 0, the real part of af is positive 



since Ej 3 < 0. This implies that when PA laser is 
tuned on resonance or on the blue side of the resonance, 
the modified two-body interaction is repulsive. On the 
other hand, when S v j < Ej hl , the real part of af is neg- 
ative. This means that when PA laser is tuned on the red 
side of the resonance by an amount exceeding \E s J i ' l * t \, 
the modified interaction becomes attractive. For the pa- 
rameters given in Table-I and assuming S v j = — 15 X KTi, 
we make an estimate of af =1 / a° l=1 — 11. From low energy 
behavior of unperturbed scattering wave functions, it fol- 
lows that Ti — k 2e+1 for I = 1. Therefore, af ~ k l as 
k — > 0, which is significantly different from the behavior 
of background scattering length aP e ~ k 21 . This clearly 
demonstrates that nonzero partial wave scattering ampli- 
tudes can be significantly modified by OFR. The modi- 
fication of p-wave scattering state can be experimentally 
observed by sending a weak probe laser beam tuned near 
J = 3 transition while keeping the intense laser beam 
(tuned near J = 1) operational as shown in Fig.l. Since 
J = 3 level can not be populated by PA transition from 
s-wave scattering state, the appearance of J = 3 line in 
PA spectra will unambiguously reveal optically induced 
p-wave Feshbach resonance. The modification can also 
be enhanced by applying another intense laser field tuned 



near J = 3 level as schematically shown with the double 
blue arrow (on the right) in Fig. 1. 

In conclusion, we have demonstrated that not only 
s-wave but also higher partial wave atom-atom interac- 
tion can be manipulated by the method of optical Fesh- 
bach resonance with an intense PA laser. We have given 
quantitative estimate of relative modification of p-wave 
scattering amplitude with a model calculation without 
hyperfine interaction. However, inclusion of hyperfine 
interaction will not alter the qualitative nature of our 
main results which are: (1) As a result of strong-coupling 
PA laser-induced large light-shifts, atoms experience dis- 
persive light force leading to modified atom-atom inter- 
action. (2) PA laser-induced modification changes the 
threshold behavior significantly. 
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